Centrality dependence of strangeness production in heavy-ion collisions as a 
geometrical effect of core-corona superposition 
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It is shown that data on strange particle production as a function of centrality in Au-Au collisions 
at \/s NN = 200 GeV can be explained with a superposition of emission from a hadron gas at full 
chemical equilibrium (core) and from nucleon-nucleon collisions at the boundary (corona) of the 
overlapping region of the two colliding nuclei. This model nicely accounts for the enhancement of 
<j} meson and strange particle production as a function of centrality observed in relativistic heavy 
ion collisions at that energy. The enhancement is mainly a geometrical effect, that is the increasing 
weight of the core with respect to corona for higher centrality, while strangeness canonical suppres- 
sion in the core seems to play a role only in very peripheral collisions. This model, if confirmed at 
lower energy, would settle the long-standing problem of strangeness under-saturation in relativistic 
heavy ion collisions, parametrized by 7s. Furthermore, it would give a unique tool to locate the 
j /' onset of deconfinement in nuclear collisions both as a function of energy and centrality if this is to 

' be associated to the onset of the formation of a fully equilibrated core. 

<N 

I. INTRODUCTION 

One of the observed features of hadron production in relativistic heavy ion collisions is the deviation from full 
chemical equilibrium of particles containing strange constituent quarks. This is described by a phenomenological 
factor 7^ s (n s being the number of strange quarks in the given hadron species) which multiplies the equilibrium 
abundance of hadrons and turns out to be generally < 1. Actually, 75 shows a mild increase as a function of centre- 
of-mass energy for central collisions from 0.65 at y/s NN ~ 4.5 to about 1 at 200 GeV [lj, according to most analyses 
0. However, even at the largest energy of Au-Au collisions, 7s turns out to be less than 1 in peripheral collisions, 
' showing a monotonically decreasing trend going from central to peripheral [H, H, 0| • 
^>0 , The idea that 75 could be the effect of superposing particle production from different sources in a single heavy ion 
collision was put forward in refs. [H, Q. In ref. [H, 0|, the multiplicities of various hadron species in Pb-Pb collisions 
at SPS were described well with this core-corona picture assuming that corona is a halo of single nucleon-nucleon 
collisions where produced particles escape the interaction region unscathed, while the core gives rise to a completely 
equilibrated hadron gas, i.e. with 7s = 1. Since strangeness production is suppressed in NN collisions with respect 
to a fully equilibrated hadron gas, while temperature is almost the same [7|, |8j), if such single NN collisions accounts 
for a significant fraction of total particle production, a global fit to one hadron-resonance gas would actually find 
7s significantly less than 1. In ref. [6j, the authors assume that a string percolation process gives rise to a large 
cluster in the core of the nuclear overlapping region and smaller clusters in the outer region (henceforth referred 
to as corona), eventually decaying into hadrons according to the statistical model ansatz. With this core-corona 
5_J ■ superposition scheme, and assuming js = 1 the authors could reproduce the centrality dependence of K/ir ratio at 
SPS and RHIC because small corona clusters suffer the so-called canonical suppression effect. 

This core-corona model has been applied to other observables. It has been found to be able to describe rapidity 
densities of charged hadrons in Au-Au collisions More detailed analysis found out that the rapidity densities 
of various hadron species as a function of centrality in Au-Au collisions at 200A GeVat RHIC as well as nuclear 
modification factors are well described with the EPOS model in a core-corona scheme[l(|. Finally, it has been taken 
into account also for analysis of J/ip production within the statistical hadronization model [ill ]. 

In this paper, we will show that this core-corona superposition nicely accounts for strangeness enhancement as 
a function of centrality observed at RHIC. The key probe which demonstrates the viability of this picture is the (f> 
meson, which, being a completely neutral particle, does not feel canonical suppression. This, as we will discuss in 
detail in Sect. IIIH favours the picture of corona as originated from NN collisions rather than small clusters hadronizing 
into a fully equilibrated hadron gas. In fact, in the statistical hadronization model, sufficiently small clusters entail 
canonical suppression for open strange particles, but not for those with hidden strangeness. 

In ref. [7J, we fitted the number of single nucleon-nucleon collisions taking place in the corona as a free parameter. 
In this paper, we will assume a definition of corona as those nucleons which undergo one collision and show that this 
succesfully accounts for the meson centrality dependence. Indications that strangeness suppression is related to the 
number of multiply colliding nucleons were found in ref. and, very recently, in ref. [l2j ; in this work, we clarify this 
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relation. 

The basic ideas and conclusions discussed here have been reported earlier in ref. [13]; in this work, we expand, 
explain and update our analysis. 

II. STATISTICAL MODEL AND CANONICAL SUPPRESSION 

Statistical model analyses [J, 0, E3] in Au-Au collisions at 200^4 GeV with mid-rapidity densities find that the 
chemical freeze-out temperature as well as the baryon chemical potential are constant throughout the accessible 
centrality range implying that the thermodynamical state of the produced matter at mid-rapidity does not depend 
on centrality at chemical freeze-out. Not so for the strangeness under-saturation parameter 7s which is found to be 
significantly less than 1 in peripheral collisions. 



It has been argued [15| that 7s < 1 is an effect of so-called canonical suppression effect. Namely, strange particles 
are suppressed with respect to their expected yield in a grand-canonical ensemble (or thermodynamic limit) because 
strangeness is exactly vanishing within a small volume, called strangeness correlation volume (SCV), which does 
not coincide with the volume of the average fireball at mid-rapidity or the global freeze-out volume of fireballs. 
Therefore, going from pp collisions to central heavy ion collisions through peripheral ones, one should observe a 
relative enhancement of strange particles due to approaching the thermodynamic limit, which is hierarchical: Q yield 
increases faster than S which increases faster than A's or kaons. Yet, although this hierarchy of enhancements has 
been observed both at SPS 16] and RHIC [13], in neither case the natural saturation expected when SCV attains 
a sufficiently large value is seen. In fact, this means that the SCV would only reach its saturation value (the one 
sufficient for the system to be essentially grand-canonical) at RHIC precisely in central collisions, where 75 ~ 1. 
Therefore, we think that canonical suppression is quite an unnatural explanation of the data, as pointed out by many 

[II [HE [II. 

There is, however, a clearcut probe to test the canonical suppression picture and this is <f> meson. It is not an open 
strange particle, thus it is not canonically suppressed, yet, being a ss state, it must be 7^ suppressed. Furthermore, 
<j> meson has almost no feeding from heavier species, so it does not suffer canonical suppression even indirectly as a 
decay product of open strange particles. It was pointed out quite early [20| that a statistical model with canonical 
suppression mechanism, i.e. with SCV as additional parameter, would have not been able to explain the deviation of 
the 4> meson yield from its grand-canonical value and this has been demonstrated in fits to NA49 multiplicities [f|. 
Recently, STAR collaboration has measured the mid-rapidity densities of <fi meson very accurately and the observed 
pattern as a function of centrality clearly shows (see fig. [T|) that these do not scale linearly with the number of 
participants, rather the ratio to pp value increases rapidly at very peripheral collisions slowly saturating thereafter. 
This non-linear increase cannot be accounted for by a variation of the chemical freeze-out temperature because this is 
constant as a function of centrality, as has been mentioned. The only way to accommodate the <f> meson behaviour in 
a statistical model fit is to introduce a 75 factor, which is found to be significantly less than 1 in peripheral collisions. 

Instead of introducing an ad hoc parameter to describe centrality dependence of <j) meson and other strange particle 
yields, we can try to explain strange phase space under-saturation as an effect superposing two different particle 
sources, as has been mentioned in the Introduction: a fully equilibrated core (i.e. a hadron gas with 75 = 1) and 
single NN collisions in the corona. The appearance of 7s in global statistical model fits is owing to the suppression of 
relative strangeness production (with respect to hadron gas in full equilibrium) in the unavoidably present single NN 
collisions in the corona. This picture would naturally account for the decrease of fitted 75 in peripheral collisions, 
where the core is smaller and the importance of single NN collisions grows. 

It should be pointed out that statistical model fits to pp collisions generally find temperatures only 10% higher than 
in heavy ion collisions at the same beam energy but consistently lower 75 pp — 0.5, i.e. about a factor 2 smaller [3]. 
This explains why, in this core-corona model, no effect is seen in global fits on centrality dependence of temperature 
but a significant dependence of 75. 

Finally, the shape of normalized <f> mid-rapidity density in fig. [T] also tells us that the corona cannot be really made 
of small clusters hadronizing into a fully equilibrated hadron gas at the same temperature of the core. Indeed, in this 
case, 75 = 1 and one could account only for the suppression of open strange hadrons but not of <j>. 

III. CORE-CORONA SUPERPOSITION 



In the following, we will introduce a very simple model which is suitable to study core-corona superposition in heavy 
ion collisions and, particularly, to estimate the relative weight, in terms of particle production, of core and corona as 
a function of centrality. 
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It has been observed [2l|,[22| that the number of charged hadrons emitted in hadron-nucleus (hA) collisions, compared 
with pp collisions at the same beam energy, scales as: 

/dN ch \ /dN ch \ Np A N P N coll + 1 

\ dy / J\ dy / pp iVPP 2 2 ' U 

where A co n is the number of collisions. Based on the above formula, the rapidity density of a given hadron species i 
in heavy ion collision at mid-rapidity is written as a sum of two contributions: 

m \ = Npc/m\ + /m 

dy J 2 \ dy / NN \ dy 



Npc/m\ /m 

2 \ dy I \ dy 



(2) 



In the above equation, Ape is the mean number of participants in the low density corona, which undergo single 
nucleon-nucleon collisions whose produced particles escape the interaction region unscathed. Since at high energy 
hadron production in neutron-neutron as well as in neutron-proton collisions closely resembles that in proton-proton 
collisions at the same beam energy, the second (approximate) equality in Eq. ((2]) holds true. The rapidity density in 
pp collisions is measured, while that in the core is assumed to be that of a completely equilibrated hadronic source. 

At very high energy (typically RHIC energies) where particle rapidity distribution is wide, the second term on the 
right hand side of Eq. ([2|) can be written Q : 

m\ -vjp) <3) 



\ d y /core \ d V 

where n, is the density of hadron species i in the fully equilibrated average fireball of freeze-out volume V and po is 
the density of fireballs per unit of rapidity at mid-rapidity. For a collision of nuclei with mass numbers A and B, the 
volume of the core at freeze-out can be estimated: 

V . m - SVo) „ , (JM - _ f (NT _ Nrc} . (4) 

In Eq. ^ Np is the number of participants, Ape is the number of participants in the corona; Vq ~ Np/2n,Q (uq being 
some initial density related to nuclear density) is the average of the initial overlap volume of the two colliding nuclei; 
SVq w Npq/2tiq is the average volume of the corona, / is a "growth factor" which takes into account the expansion of 
the system between the initial overlap time and freeze-out. It should be stressed that the essential point here is the 
proportionality between the core volume and the number of participants of the core obtained as difference between 
the total number and that of the corona; indeed, as we will see, the knowledge of the parameters / and no is not 
needed for our purpose. 
Plugging Eq. @ into © we obtain: 



(N P - Np C ) ( ^1 ) . (5) 



Since charged particle multiplicity scales linearly with the number of participants in heavy ion collisions, we are led 
to conclude that the factor fpo/2riQ in Eq. ([S]) is independent of A^p, i.e. independent of centrality, provided that 
(drii/dy} COIC is in turn independent of centrality. As we will see, this condition holds true because of the observed 
independence of freeze-out parameters on centrality. 

Putting Eq. (J5J) into @ and dividing both sides by Ap (dn/dy) /2, we obtain: 



pp ' 

dNi \ /dm' 



R _ Z \^/aa _ 2/po W core ( A PC \ A PC _ A PC / A PC \ 
A Np /m\ 2n /d^\ V A P y A P A P ^ N P J 



dy / \ dy 



(6) 



The factor A embodies all unknown parameters and depends linearly on particle density in the core. If this factor 
was independent of Ap , the whole centrality dependence of the Ra ratio would be given by quantity related to the 
geometry of the collision. Indeed, we now need to define the NN collisions which form the corona, or, in other words 
to define Ape- 
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FIG. 1: Ratio Ra (see text for definition) for the <f> meson in Au-Au collisions at = 200 GeV. The data points are from 

STAR 23] while the lines are calculated according to the Eq. ((8| by fixing A in the 2nd most central and in the mid-peripheral 
bin. 



Based on Eq. for hA collisions, where the impinging and target nucleons collide only once, we choose to define 
A^pc as the number of nucleons colliding only once. Explicitly: 



N PC = Ng + Ni 



(7) 



in which N s ' are the numbers of such singly colliding nucleons from nuclei A and B respectively 1 . 

Armed with the definition ([7]) of single collisions in the corona, and with the equation (JSJ), we are now in a position 
to test the model on the data by calculating Np and Npc with the Glauber model. 



IV. DATA ANALYSIS AND RESULTS 

We have emphasized that the <f> meson provides us with an excellent probe to test models of strangeness suppression, 
and so does it for the core-corona model and the Eq. ([6]), that we rewrite here: 



Ra 



dNj 
dy 



AA 



dNj 
dy 



A + ^£(i-A). 
N P V ' 



(8) 



pp 



Since <fi is immune from finite size effects, i.e. canonical suppression, the factor A is independent of centrality because 
chemical freeze-out temperature is found to be centrality- independent. We can then fix the factor A from the data in 
one centrality bin and see how well the normalized yield Ra is described at other centralities by calculating Np and 
A^pc with a Glauber Monte-Carlo simulation. 

For our Glauber Monte-Carlo calculation we have implemented essentially the same algorithm used in ref. [24j. Our 
calculated number of participants in the corona according to formula is shown in Fig. O In Fig. [T] Ra for the <f> 
meson measured at RHIC at V^TVAf = 200 GeV [23| is compared with the theoretical calculation. The two lines are 



1 In ref. [ 1 3t | we have used a different definition of -/Vpc, i.e. twice the minimum between Ng and Ng . We find that the present definition 
is conceptually more satisfactory and better motivated experimentally. However, they eventually lead to very similar results. 
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calculated according to the Eq. ([8]) by fixing the factor A from two different centrality bins. We can see that in both 
cases the agreement between model and experiment is excellent. 

40.0 i , , . , , , . 1 
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FIG. 2: Number of nucleons Ape undergoing one collisions denning the corona as a function of Np in Au-Au collisions at 
yG NN = 200 GeV. 

On the left panel of Fig. [3J similar curves are shown for A, A, E ± and Q + f2 compared with Ra's measured at 
RHIC [17] at the same beam energy. The factor A is fixed to the 2nd most central bin 2 and, as expected, the model 
overshoots the data in the most peripheral collisions. This is because the factor A is not in fact independent of 
centrality and should decrease in very peripheral bins due to strangeness canonical suppression. We note in passing 
that the difference between the curve and the points gives then quantitative information about the volume of the 
average fireball at mid-rapidity. 

In Fig. 3] we show the fraction of particle production from core and corona as a function of centrality for different 
particle species. This plot shows that the hierarchy of enhancement slopes observed at RHIC is nicely reproduced in 
this approach. According to Eq. (JSJ), the higher is A, the steeper is the increase of Ra as a function of Np. Since A 
depends on particle species through the ratio: 

\ u I core 

the observed hierarchy simply reflects the hierarchy of ratios of mid-rapidity densities in heavy ion to NN collisions, 
therefore Q > S~ > S + > > A > A as it results from Fig. [4] 

The model can be further tested with the RHIC data directly with the Eq. (0 in the statistical hadronization model 
framework (see ref. The hadron radiation from the coro na p art can be estimated by taking the experimental 

rapidity densities of different hadron species i in pp collisions [23L [26l . |27| measured by the STAR collaboration and 
multiplying these with the number of corona participants Npq. This number can be determined in a twofold way: by 
fitting it as a free parameter or calculating it directly from the Glauber model as the total number of singly colliding 
nucleons as in Eq. 0. In the first case we will have 4 free parameters that must be fitted to the measured rapidity 
densities at different centralities while in the second case we have 3 free parameters only. The fitted and calculated 
Npc as a function of Np at different centralities are shown in fig. [5l One can see that the fitted and calculated 
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FIG. 3: LEFT: Ratio Ra (see text for definition) for hyperons measured in Au-Au collisions at \/*nn = ^00 GeV. The data 
points are from STAR [ItJ while the lines are calculated according to the Eq. (|8]) by fixing A in the most central bin. 
RIGHT: The same quantity Ra for tv + , K + and p in Au-Au collisions at v^nn = 200 GeV. The data points are from STAR 
[25l ] while the lines are calculated according to the Eq. (|8} by fixing A to the 2nd most central bin. 




FIG. 4: Fraction of produced particles coming from core (upper lines) and corona (lower lines) as a function of centrality. 



number of single collisions agree very well with each other in peripheral and semi-central collisions while in central 
collisions the relative emission from the corona is too small for us to reliably fit the number of single collisions. We 
also remark that hadron yields are described better in this two-component formalism than the conventional statistical 
hadronization model analysis [l[ with the parameter. 
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FIG. 5: Number of corona participants iVpc at different centralities. The square dots denote the values calculated from Glauber 
model while the round ones denote the values arising from fitting iVpc as a free parameter. The round symbols are shifted 5 
units of iVp rightward for clarity. 



V. SUMMARY AND CONCLUSIONS 

In summary, we confirm our early finding [l3| that strangeness enhancement from peripheral to central relativistic 
heavy ion collisions at RHIC at \^s NN = 200 GeV can be well described by a model where particle production arises 
from two sources: a fully equilibrated core at a temperature of ~ 160 MeV as in the statistical model and an outer 
region of single NN collisions, called corona. Since in NN collisions relative strangeness production is suppressed with 
respect to a fully equilibrated grand-canonical core, the observed enhancement going from peripheral to central stems 
from the increased weight of the core with respect to corona. The enhancement is thus mainly a geometrical effect 
and it is not driven by the conservation of net strangeness within small regions in the core itself (canonical suppression 
picture): this effect shows up only in very peripheral collisions where the whole core's volume is presumably small. 

The key probe for our argument is the (f> meson which is immune from canonical suppression and has essentially no 
feeding from high-lying resonances. The observed rise of relative 4> yield as a function of centrality is an unambiguous 
signal that the increased production of strangeness from peripheral to central collisions is not an effect of strangeness 
conservation nor can it be explained by an increase of temperature, which is found to be constant throughout. The 
4> enhancement also favors the idea that corona consists of independent NN collisions rather than small clusters 
hadronizing in full chemical equilibrium. In the latter case, with a hadronization temperature at the same value of 
around 165 MeV (as confirmed by analysis of pp collisions), <f> normalized production would be flat as a function of 
centrality. 

The success of this description indicates, as pointed out, that a fully equilibrated core is formed at any centrality 
in Au-Au collisions at y/s NN — 200 GeV. Arguably, this should occur also at lower energies and so core-corona 
superposition is likely to explain strangeness enhancement observed by NA57 experiment in Pb-Pb collisions at 
V^jvat = 17.3 GeV. It becomes then crucial, in order to locate the onset of full chemical equilibrium in the core, to 
study the production of hyperons and chiefly <f> mesons as a function of centrality in the energy range from few GeV's 
to 20 GeV. When, at some low energy or centrality, the data will start to be well reproduced by a fully equilibrated 
core, that might be be the point where deconfincment has occurred. 

Finally, the superposition of core and corona appears to be a general feature of relativistic heavy ion collisions which 
should be taken into account for the analysis of all observables besides particle chemistry. This has been pointed out 
in refs. [jjl [To| and in a very recent analysis (28j. The definition of the corona is not unique and our results indicate 
that it should be based on NN collisions rather than local density [l(| El • In this regard, a global fit to particle 
spectra [H| also favors our definition of corona as being made of those nucleons undergoing one collision. 
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